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ABSTRACT 

We present an analysis of a new, detached, double-lined eclipsing binary system with 
K7 Ve components, discovered as part of the University of New South Wales Extrasolar 
Planet Search. The object is significant in that only 6 other binary systems are known 
with comparable or lower mass. Such systems offer important tests of mass-radius 
theoretical models. Follow-up photometry and spectroscopy were obtained with the 
40-inch and 2.3m telescopes at SSO respectively. An estimate of the radial velocity 
amplitude from spectral absorption features, combined with the orbital inclination 
(83.5°) estimated from lightcurve fitting, yielded a total mass of M tota / = 1.041 ± 0.06 
M and component masses of M A = 0.529 ± 0.035 M and M B = 0.512 ± 0.035 M . 
The radial velocity amplitude estimated from absorption features (167±3 km s _1 ) was 
found to be less than the estimate from the H a emission lines (175±1.5 km s _1 ). The 
lightcurve fit produced radii of R A = 0.641 ± 0.05 R and R B = 0.608 ± 0.06 R , 
and a temperature ratio of Tb/Ta — 0.980 ± 0.015. The apparent magnitude of the 
binary was estimated to be V — 13.9 ±0.2. Combined with the spectral type, this gave 
the distance to the binary as 169 ± 14 pc. The timing of the secondary eclipse gave a 
lower limit on the eccentricity of the binary system of e ^ 0.0025 ± 0.0005. This is the 
most statistically significant non-zero eccentricity found for such a system, possibly 
suggesting the presence of a third companion. 
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late-type - stars: individual UNSW-TR-2 



1 INTRODUCTION 

Low-mass stars (M < Mq) make up more than 70% of 
the stellar population, but their intrinsic faintness means 
that relatively few have been studied in detail. In recent 
years, much work has been don e on the theoretical mod - 
elling of such stars (see review bv lChabrier fc Bara ffc 2000), 
but progress has been hindered by a lack of observed low- 
mass stars with sufficiently accurate parameters. In order 
to distinguish between the various models, uncertain ties of 
less than 2 -3% in the masses and radii are required (| Ribasl 
I2005l . l2003h . This can most easily be achieved with detached, 
double-lined eclipsing binary stars, with best results c oming 
from systems with similar components i|Ribasll2005l) . Cur- 
rently, only six such systems with masses of 0.6 Mq or be- 
low ar e known. A number of recen t studies of these sys- 
tems (lLopez- Morales fc Ribasl 120051: iTorres fc Ribasl liooj 
iMace roni fc Montalb an 2004; Ribas 2003J) have revealed dis- 
crepancies between the observations and the latest theoreti- 
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cal models, which in particular appear to underestimate the 
radii by as much as 10-15%. This is possibly due to the en- 
hanced magnetic activity produce d by the hig h rotational 
velocities of stars in close binaries (|Ribasll2005l ). There is a 
clear need for more of these systems to be discovered and 
studied, so that the physics of low-mass stars may finally be 
understood. 



In this paper we report the discovery and preliminary 
analysis of a new detached, double-lined eclipsing binary 
system with K7 Ve components. The system was detected 
in data obtained for the Universit y of New South Wales 
(UNSW) Extrasolar Planet Search ((llidas et al.ll2005l U The 
eclipses are grazing, and their depth was diluted by light 
from a brighter, blended star. The resulting shallow eclipses 
led to this object's selection as as a transiting planet candi- 
date, its true nature only being revealed by follow-up obser- 
vations. 
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Table 1. Infrared magnitudes of UNSW-TR-2 and its close neigh- 
bour from the 2MASS catalogue. 





TR-2 


Neighbour 


J 


10.9 


10.7 


H 


10.2 


10.2 


K 


10.1 


10.1 



2 OBSERVATIONS 
2.1 Photometry 

The data leading to the detection of UNSW-TR-2 
(«J2000 = 18 h 30' 52'.'3, S J20 oo = +7° 9' 27") were obtained 
using the 0.5 m Automated Patrol Telescope (APT) at Sid- 
ing Spring Observatory (SSO), NSW, Australia. The object 
was found in a field centred on the open cluster NGC 6633. 
During three seasons (2002-2004) we observed this field on 
79 nights (about 3500 images). For most of the present anal- 
ysis, we have used only the best quality data (highest ca- 
dence and lowest photometric noise), consisting of 997 im- 
ages from 11 nights in July and August 2002. On these nights 
we observed through a Johnson V filter at a cadence of about 
16 per hour. The remaining data on this field have one quar- 
ter the cadence (we were observing three other fields in par- 
allel) and are affected by significant systematic trends on 
many nights. For the period determination (sec. 13. ip . we 
used the entire dataset. Lightcurves were extracted using 
a custom-built automated reduction pipeline. For more in- 
formation a bout the tel e scope , the data and the reduction 
process, see iHidas et al.1 (|2005l ). 

Initial inspection of the lightcurve of UNSW-TR-2 re- 
vealed eclipses with a depth of ~25 mmag, a duration of 
2 hours, and a binary period of ~ 2.1 days. The out-of- 
eclipse magnitude is V — 12.0 ± 0.1, but inspection of a 
Digitised Sky Survey (DSS j3 image revealed that the APT 
photometry aperture contains at least two blended stars of 
similar brightness. The APT has 9.4-arcsecond pixels, and 
in a crowded field such as this one (at galactic latitude +8°), 
such blending is common. The photometry aperture is 3 pix- 
els in radius. 

In June 2004 we obtained images of UNSW-TR-2 at 
higher spatial resolution in V and I using the 40-inch tele- 
scope at SSO. The images were sampled by 0'.'6 pixels. We 
observed one complete eclipse (in ~ 3" seeing) , with a depth 
of 0.14 mag. The eclipse has the same depth and shape 
(within the photometric errors) in the two colour bands, 
suggesting that the two stars are of similar temperature. 

The binary system is in fact the second-brightest object 
in the APT photometry aperture. It is 1.48±0.03 mag fainter 
than the brightest o ne, and separated from it by ~10". The 
Tycho-2 catalogue l|H0g et al l koOOh gives a (Johnson) V 
magnitude of 12.4 ±0.2 for this bright neighbour. The com- 
bined apparent magnitude for the binary system is therefore 
V = 13.9 ±0.2. 

The Digitised Sky Survey was produced at the Space Telescope 
Science Institute under U.S. Government grant NAG W-2166, 
based on photographic data obtained using the Oschin Schmidt 
and UK Schmidt telescopes. 



2.2 Spectroscopy 

2.2.1 Spectra Obtained 

Spectra were obtained using the Double-Beam Spectrograph 
(DBS) on the 2.3 m telescope at SSO. During our initial 
follow-up run in June 2004, we obtained two spectra of 
UNSW-TR-2 at the same epoch (Figure [TJ). These covered 
the range 5800-6600 A with 60 km s" 1 resolution (FWHM, 
as measured from the sky lines in the spectra). 

Visual comparison of the flux-calibrated spectra with 
templates from the UVILIB spectral library (| Pickles! [l998h 
suggested that both stars are of type K7 V. With addi- 
tional spectra reaching further into the red, our estimate of 
the spectra l type could be refined using the TiO absorption 
band s (e.g. iKenvon fc Fernandez-Castrolll987l ; ISchild et al.l 
1992). One TiO band is present in our existing data, but the 
absorption is very weak. 

All of the absorption lines have two components, the 
smaller of which has 60-80% of the depth of the larger com- 
ponent, pointing to two similar stars. Prominent H a emis- 
sion is also present, with the two components at a similar 
strength-ratio to the absorption-line pairs (Fig. [TJ . 

The expected absolute magnitude of the two stars com- 
bined, bas ed on their spectral type (K7 V), is Mv = 7.6 
{Coxl 12000). Combining this with our estimate of the sys- 
tem's apparent magnitude (V — 13.9 ± 0.2), and assuming 

1 mag kpc -1 of extinction, puts the binary at a distance of 
169 ± 14 pc. The template best matching the spectrum of 
the bright neighbouring star is that of a G8III giant, with 
Mv = 0.8. Assuming the same extinction gives a distance 
of 1200 ± 70 pc. Therefore this neighbour is not physically 
associated with the binary system. 

2.2.2 Radial Velocity Estimation 

We fitted for the two stellar components simultaneously us- 
ing VPFITl). For both spectra, five absorption-line pairs 
with high signal-to-noise were individually fitted with 
blended Voigt profiles. When given the instrumental reso- 
lution, VPFIT fits for the central velocity, column density 
and Doppler broadening of each line and returns associated 
uncertainties based on the covariance matrix at the best- 
fit. Each absorption-line pair is surrounded in our spectra 
by weak absorption features, some of which are partially 
blended with the pair in question. To allow for this, we in- 
serted extra absorption lines into VPFIT as appropriate, 
again allowing all line-parameters to vary. We found that in 
order to obtain reasonable fits it was also necessary to allow 
for the presence of a final weak subsidiary line internal to 
the main absorption-line pair. In all cases, this procedure 
produced a reasonable fit to the data (x 2 ~ 1)- 

A small adaptation of VPFIT allows for the fitting of 
emission lines with blended Gaussian profiles. In this way 
we fitted the two H a emission-line-pairs. No subsidiary lines 
were required to obtain a good fit. 

Taking orbital phase into account, and assuming a cir- 
cular orbit, the total velocity-amplitude (Ka ± -Kb) sin i is 
167±3 km s _1 when measured from the absorption features, 
and 175 ± 1.5 km s _1 when measured from the H a emission 

2 http:/ /www. ast.cam.ac.uk/~rfc/vpfit. html 
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Figure 1. The spectrum of UNSW-TR-2, obtained at phase 0.213 using the 2.3 m telescope at SSO. The inset shows a close-up of the 
H Q lines. 



lines, where i is the orbital inclination of the system (sub- 
sequently determined by lightcurve fitting, as described in 
section I3.2|l . An apparent disagreement between velocities 
derived by these two methods has been noted in other sys- 
tems (M. Lopez- Morales, private communication). The ab- 
sorption velocity is considered to be more reliable due to un- 
certainties as to the spatial distribution of the source of H Q 
emission relative to the photometric centre of each star. We 
therefore adopt the value of (Ka+Kb) sin i — 167±3 km s _1 
for the velocity amplitude of our system. We were unable to 
improve upon this estimate using the two-dimension a l cros s- 
correlation software TODCOR (|Mazeh fc Zuckerl (|l994l )) 
with a variety of different template spectra. This is possi- 
bly due to the relatively low resolution of our spectra when 
compared to those of similar objects for which TODCOR has 
been successfully implemented (e.g. lLopez-Morales fc Ribasl 
l2005h . 

In March 2006 we obtained an additional spectrum with 
the same instrument and at the same resolution. This one 
was taken close to an eclipse, therefore the lines are com- 
pletely blended and the separation of the absorption lines 
cannot be reliably measured. However, comparing the cen- 
tral wavelength of each blended line with the wavelengths 
of the separated components in the original spectra allows 
us to constrain the difference between the two masses. This 
was done for the five high-signal-to-noise absorption lines 
in both of the original spectra, leading to a constraint of 
AM = Ma - Mb = (0.016 ± 0.03)Mt otai . 



3 ANALYSIS 



3.1 Period and Eccentricity 

Phase-folding the entire APT lightcurve (with observed 
eclipses up to two years apart) allows a precise determi- 
nation of the orbital period. The period giving the smallest 
scatter in the folded lightcurve is 2.11674 ± 0.00002 days. 

The epoch of the primary eclipse minimum (to) was de- 
termined from the 40-inch lightcurve. The secondary eclipse 
occurs at a phase of 0.5016±0.0003 (measured from the APT 
lightcurve), 5± 1 minutes later than expected for a perfectly 
circular orbit. Assuming that the major axis of the orbit lies 
in the plane of the sky gives a lo wer limit on the orbit's 
eccentricity of e > 0.0025 ± 0.0005. IZahn fc Bouchetl ||l989h 
predict that a binary consisting of two 0.5 Mq stars with a 
period of about 7 days and an initial eccentricity of 0.3 will 
have been circularised by tidal interactions to e ~ 0.005 by 
an age of 10 6 years, and will have reached e ~ 0.004 after 
10 10 years. However, the shorter period of our system should 
lead to increased circularisation and lower e. Although our 
lower limit may be compatible with a simple, main sequence 
binary system, it suggests the possibility that our system ei- 
ther is very young, or more likely, has at least one other 
companion affecting the timing of the eclipses. This is the 
most significant evidence for non-zero eccentricity in a low- 
mass eclipsing binary system detected to date. However, we 
continue to assume e = for this work. 
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Figure 2. The four eclipses used in the parameter estimation. Top: The primary eclipse observed in both V and / with the ANU 40-inch 
telescope. Bottom: Both eclipses from the APT lightcurve, phase folded and binned at 5-minute intervals. The solid line shows the best-fit 
lightcurve. 



3.2 Lightcurve Fitting 



We fitted models (based on equations by ISackettJ (| 1999T ) . 
modified to apply to binary stars) to the lightcurves using 
a chi-squared minimisation algorithm to estimate properties 
of the system. The APT lightcurve was folded and binned at 
5-minute intervals. All the data (one high-resolution primary 
eclipse in each of V-band and /-band from the 40-inch, and 
both eclipses in V from the APT) were fitted simultaneously. 
The accuracy of our results is limited by the lack of high- 
spatial-resolution observations of the secondary eclipse, and 
also by the presence of the nearby brighter star in the APT 
aperture. The epochs of the primary and secondary minima, 
the period, and the velocity amplitude were fixed at their 
previously measured values whilst allowing the orbital incli- 
nation (i), the radii of both stars (Ra and Rb), and their 
temperature ratio (Tb/Ta) to vary. Also required is the 
fraction F of the total flux in the APT aperture produced 
by our object, estimated as 0.19 ± 0.02 from the 40-inch im- 
ages and allowed to vary within this error-range. Finally, we 
assume solar-metallicity limb-dar kening coeffic ients of 0.9 
and 0.7 in V and I respectively (|Claretl H998T I. The best- 



fit lightcurve is shown in Figure [2] and the parameters are 
summarised in Table [2] 

3.3 Error Estimation 

To estimate the uncertainties in our lightcurve fit we plotted 
Ra vs. Rb reduced-chi-squared contour plots for ranges of i, 
Tb/Ta and F either-side of the best-fit parameters (Fig. [3}. 
The contours on each Ra vs. Rb plot are ellipses, tightly 
constraining the sum of the radii but not their difference. 
Altering the other parameters from their best-fit values de- 
creases the size of the contours whilst moving their centre. 
Increasing i leads to a higher Ra but lower Rb, whilst in- 
creasing Tb/Ta has the opposite effect: lower Ra and higher 
Rb- We found that the flux ratio F is constrained relatively 
tightly and therefore contributes little to the uncertainties 
of the other parameters. Having scaled the errors such that 
the best-fit reduced-chi-squared has a value of one, the la 
errors for the remaining four parameters were estimated by 
the range of the appropriate chi-squared contour over the full 
4-dimensional error ellipsoid. Further (minor) contributions 
to these errors come from the uncertainties in the velocity 
amplitude and limb-darkening coefficients. Because the es- 
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Figure 3. The Ra vs. Hb chi-squared contour plot around the 
best-fit solution to the lightcurve. The other parameters are set 
to their best-fit values. The three contours represent the la, 2a 
and 3a confidence regions, but to estimate the true uncertainties 
Tg /Tj4 and i must also be varied. 
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Figure 4. Mass vs. radius for all known components of double- 
lined eclipsing binary stars with M < 0.8 Mq. UNSW-TR-2 is 
shown as triangles. The solid line is an ap proximation to a 300 - 
Myr isochrone calculated with models from l Baraffe etaT. ] lll99St) . 
The figure is based on Fig. 1 in lRibasI ^OOcfT 



timates of the two radii are not independent, we also used 
the contour plots to derive an estimate of AR = Ka — Kb- 



4 RESULTS AND DISCUSSION 

The parameters we have derived for UNSW-TR-2 from both 
the lightcurve fitting and spectral analysis are shown in Ta- 
ble [2] The derived masses are consistent with our earlier 
estimate (based on spectral shape) of two K7 V stars. The 
star estimated to have a higher mass, radius and tempera- 
ture corresponds to the component identified in section [2.2.1l 
as having stronger absorption and emission features. 

In Figure [3] we plot the masses and radii of all the 
known components of detached, double-line eclipsing bina- 
ries with masses below 0.8 Mq, including UNSW-TR-2. We 



Table 2. Physical properties of UNSW-TR-2, derived from 
radial-velocity measurements (sec. 12. 2.21 1 and lightcurve fitting 
(sec. l3T2t. 



Period (days) 
to (HJD) 
Separation (AU) 
Inclination (degrees) 
M total (Mq) 



2.11674 ± 0.00002 
2453171.0371 ± 0.0005 
0.0327 ± 0.0007 
83.49 ± 0.2 
1.041 ± 0.06 



M A (M q ) 
M s (M ) 
AM (M Q ) 
Ra (Rq) 
Rb (R©) 

AR (Rq) 

T b /Ta 



0.529 ± 0.035 
0.512 ± 0.035 
0.017 ± 0.03 
0.641 ± 0.05 
0.608 ± 0.06 
0.033 ± 0.06 
0.980 ± 0.015 



also plot a theoretic al main-sequen ce isochrone d erived by 
iBaraffe et all §998). As noted by iRibasI (|2005l ) and ref- 
erences therein, stars with masses greater than ~0.3 Mq 
seem to fall systematically above the theoretical line, with 
their radii being underestimated by ~10%. The two stars 
in UNSW-TR-2 appear to agree with this trend, although 
a more precise determination of their masses and radii is 
clearly required. 

It is likely that one or both of our stars lies in the 
previously unpopulated region between 0.5 and 0.6 Mq, 
which increases their importance in constraining the mass- 
radius relation of low-mass stars. The apparent magni- 
tude of our system (V ~ 13.9) is similar to that of the 
0.6 Mq eclipsing binary GU Bootis (V ~ 13.7), for which 
lLopez-Morales fc Ribasi (|2005T ) derived the masses and radii 
to within the 2-3% accuracy required for a detailed com- 
parison with theory. With echelle spectroscopy and further 
photometry, the properties of UNSW-TR-2 could therefore 
be similarly constrained. 



4.1 UNSW-TR-2 as a target for a planet search 

Eclipsing binary systems with dwarf components are inter- 
esting targets for photometric planet searches. One main 
advantage is that the orbits of planets in such a system are 
likely to be near-coplanar with the binary, making transits 
more likely. Additionally, non-transiting planets may also 
be detected via the variations they induce in the ti ming of 
the binary eclipse minima (e.g. lDovle fc Deedl20o4 ). These 
motivate d an extensive sea rch around the M-dwarf binary 
CM Dra jPovle et alJl200Ch . 

In the case of UNSW-TR-2, transits of a Jupiter-sized 
planet in front of one of the stars (when they are not eclips- 
ing) would be 1.3% deep. A planet in a co-planar orbit 
about the binary (i = 83.5°) would have to be within 
r max ~ 0.04 AU of the centre of mass to transit. This is 
only 2.5 times the binary's semi-majo r axis, but according 
to simulations by iHolman fc Wiegertl (|l999l ). even such a 
small orbit may be stable around an equal-mass binary with 
a circular orbit. Planets with more favourable orbital incli- 
nations (relative to the line of sight) may also transit. For 
example, i p = 85° gives r max = 0.05 AU, and i p = 87° gives 
r max = 0.09 AU. 
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The approximate range of eclipse-timing variations is 
r = qr/c, where q is the planet-to-binary mass ratio, r 
is the planet's orbit al radius, and c is the speed of light 
|Dovle fc Deed [2004'). For a Jupiter-mass planet (q ~ 10~ 3 ) 
in a 1 AU orbit, this gives r ~ 0.5 sec. The detection of plan- 
ets via this method would therefore require second-precision 
timing of eclipses over timescales of years. 



5 CONCLUSIONS 

We present a preliminary analysis of a new double-lined 
eclipsing binary system with K7 Ve components. The system 
was discovered as part of an ongoing transit search for extra- 
solar planets. The system is important in that it is only the 
seventh low-mass binary known, and because it allows pre- 
cise determinations of the masses and radii of these low-mass 
stars. Using the initial photometric data, and follow-up pho- 
tometry and spectroscopy, we derive the physical parameters 
of the system. Our results are consistent with previous work, 
confirming a departu re in the mas s -radiu s plane from the 
theoretical models of iBaraffe et al.l (|l998l ). Enhanced mag- 
netic activity — due to the high rotational velocities of stars 
in close bina ries — may be responsible for this discrepancy 
l|Ribasll2005h . 

We find statistically significant evidence — in fact the 
most significant in any known low-mass eclipsing binary sys- 
tem — for a small but non-zero eccentricity in this new 
binary system, suggesting the possibility of an additional 
companion. 
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